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Applying electric fields is an attractive way to control and manipulate single particles or molecules, 
e.g., in lab-on-a-chip devices. However, the response of nanosize objects in electrolyte solution to 
external fields is far from trivial. It is the result of a variety of dynamical processes taking place 
in the ion cloud surrounding charged particles and in the bulk electrolyte, and it is governed by 
an intricate interplay of electrostatic and hydrodynamic interactions. Already systems composed of 
one single particle in electrolyte solution exhibit a complex dynamical behaviour. In this review, 
we discuss recent coarse-grained simulations that have been performed to obtain a molecular-level 
understanding of the dynamic and dielectric response of single particles and single macromolecules 
to external electric fields. We address both the response of charged particles to constant fields 
(DC fields), which can be characterized by an electrophoretic mobility, and the dielectric response 
of both uncharged and charged particles to alternating fields (AC fields), which is described by a 
complex polarizability. Furthermore, we give a brief survey of simulation algorithms and highlight 
some recent developments. 


I. INTRODUCTION 

Dispersions of nanoparticles in electrolyte fluids are 
ubiquitous in everyday life. Prominent examples are pro¬ 
teins and DNA in aqueous environment. The study of 
such systems is not only of technological interest for the 
development of advanced materials, but also of funda¬ 
mental interest for our understanding of life science and 
biophysics. These materials have a large inherent com¬ 
plexity. Already the simplest system contains at least 
three components: large macromolecules or particles (so¬ 
lutes), small electrolyte ions, and solvent molecules. The 
interactions between the components include short-range 
interactions, such as excluded-volume and van der Waals 
interactions, and long-range interactions, such as the 
electrostatic and hydrodynamic interactions. The in¬ 
terplay of these interactions determines the equilibrium 
and dynamic properties of the system. Due to its multi- 
component nature, a wealth of parameters can be used 
to control material properties, e.g., the surface charge 
density, the salt concentration, the dielectric constant of 
the solvent. External perturbations can be used to ma¬ 
nipulate the behaviour of the system. Since charges are 
involved, electric fields are particularly efficient. 

In this review, we will focus on two important systems: 
Colloidal dispersions and polyelectrolyte solutions. Col¬ 
loidal particles are solid objects with sizes ranging from a 
few nanometers to micrometers Hi- In electrolyte solu¬ 
tion, an electric double layer forms near the solid/liquid 
interface and plays an important role in determining the 
dynamics. Polyelectrolytes are charged polymers, with 
DNA and protein being the most prominent examples 
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3. 4], In comparison to colloids, which have no internal 
degree of freedom, polyelectrolyte chains can change their 
conformation from a coil to a compact globule, and to ex¬ 
tended rodlike structures. Counterions can bind to the 
polyelectrolyte backbone in the case of strongly charged 
chains, or form a diffusive layer. In both systems, it is 
important to consider the charged particle and its sur¬ 
rounding electric double layer as a whole. 

Computer simulations have become a widely accepted 
approach to studying complex systems, complementing 
the well-established theoretical and experimental meth¬ 
ods. For charged particles in electrolyte solutions, the 
main difficulty lies in the multi-scale nature of the sys¬ 
tem. Taking colloidal particles as an example, the length 
scales range from the size of the water molecules (around 
10” 10 nr), to the size of the colloidal particle in the or¬ 
der of 10~ 6 m. Likewise, the dynamics of the system also 
involves processes with characteristic time scales span¬ 
ning several orders of magnitude. For example, the char¬ 
acteristic diffusion time (the time to diffuse along one 
molecule/particle diameter) for small salt ions and large 
colloids is of the order of 10 picoseconds and seconds, re¬ 
spectively. In principle, one can use molecular dynamics 
simulation with atomistic detail Q, but even with cur¬ 
rent computer resources, the accessible time and length 
scales are still limited. One possible solution is to use 
coarse-grained simulations, which allow one to access 
larger length and time scales at the expense of losing 
fine scale details. 

In this paper, we review recent coarse-grained simula¬ 
tion studies of charged colloids and polyelectrolyte chains 
in electrolyte solution. We mostly focus on the case of 
one single particle and its response to a weak and spa¬ 
tially homogeneous, albeit possibly time-dependent ex¬ 
ternal electric field (linear response regime). We dis¬ 
cuss the physical mechanism behind various dynamic pro- 
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cesses driven by the external electric fields. Since there is 
a vast amount of theoretical and experimental literature 
on this subject, we shall not be comprehensive, but we 
will focus on presenting simulations that illustrate im¬ 
portant physical mechanisms. We apologize in advance 
that our reference list is far from complete. 

This review is organized as follows: We start with a 
brief discussion of simulation techniques in Section E! 
Then we turn to reviewing simulation studies of particles 
and polyelectrolyte chains under constant electric fields 
(DC fields) and alternating electric fields (AC fields) in 
Section EH and hyi respectively. We conclude in Sec¬ 
tion lYlwith a brief summary and perspective. 


II. SIMULATION METHODS 

When studying the dynamics of a system that includes 
large particles (colloids or polyelectrolytes) and small sol- 
vents/ions, one needs to consider both the electrostatic 
and hydrodynamic interactions. From a simulation point 
of view, modeling such a system is a challenging task, be¬ 
cause both the electrostatic and hydrodynamic interac¬ 
tions are long-range. Taken separately, each interaction 
has been well studied in the literature [M3- Methods 
can be grouped into the categories of “implicit” and “ex¬ 
plicit” methods, SOS, depending on whether the small 
components (solvents and ions) are simulated explicitly 
or replaced by effective interactions between large com¬ 
ponents. In this review, we shall focus on explicit meth¬ 
ods. The inclusion of small ions and solvent particles re¬ 
quires more computational resources, but in many cases, 
this sacrifice is necessary for a proper treatment of the 
dynamics. For example, implicit methods that treat hy¬ 
drodynamic interactions, e.g., at the level of an Oseen 
mobility matrix, are not capable of accounting for the 
effect of finite Reynolds numbers, and difficult to com¬ 
bine with complex boundary conditions. On the other 
hand, implicit methods that replace the effect of small 
ions by screened electrostatic potentials do not capture 
the multitude of complex dynamical processes involved 
in the electrophoretic or dielectrophoretic response of a 
particle to an electric field. In the following, we first men¬ 
tion some classical approaches which are commonly used 
in mesoscopic simulations, and then highlight a few new 
developments. 

One popular method for treating Coulomb interactions 
is Ewald summation, which applies to point charges with 
periodic boundary condition [12| ■ The idea is to split the 
interaction into two contributions: one is short-ranged 
and has a cutoff, the other one is a smooth function which 
can be calculated efficiently in Fourier space. After opti¬ 
mization, Ewald summation scales as 0(J\f 3 ^ 2 ), where J\f 
is the total number of point charges. Several methods im¬ 
prove the scaling to 0(AflnAf) by computing the Fourier 
part on a grid using fast Fourier transform. Notable ex¬ 
amples are Particle-Particle-Particle Mesh (P3M) E3- 
fl5j |. Particle Mesh Ewald (PME) [16], and Smooth Par¬ 


ticle Mesh Ewald (SPME) E3- Methods that scale lin¬ 
early 0(Af) with the system size have also been proposed 
in the literature, such as the Fast Multipole Method 
ITsI , and methods based on the Maggs approach 01)1421] . 
These methods require an expensive computational over¬ 
head or have a large prefactor in the linear scaling; thus 
Ewald-based methods are still the most common choice 
in mesoscopic simulations. Some of the methods men¬ 
tioned here have been combined into a parallel library 
ScaFaCoS [22j, HH , which is freely available. 

Mesoscale methods for simulating fluid dynamics are 
based on one simple observation: As long as one is mostly 
interested in hydrodynamic effects, the solvent dynam¬ 
ics can be replaced by an artificial dynamics as long as 
the relevant conservation laws (mass, charge, momentum, 
etc.) are satisfied. Therefore, one can design simple fluid 
models that can be simulated at low computational cost. 
Popular examples in the literature include the Lattice 
Boltzmann (LB) method, Dissipative Particle Dynamics 
(DPD), and Multi-Particle Collision Dynamics (MPCD). 
DPD is a particle-based method that uses soft potentials 
and pair-wise interactions [HI, HH . LB is a lattice-based 
method which solves a linearized Boltzmann equation in 
a fully discretized fashion I n MPCD, the interac¬ 

tions between particles are performed by sor ting particles 
in cells and followed by local operations [28 - l30j . Besides 
these mesoscopic simulation methods, one can also solve 
the Navier-Stokes equations numerically, for example, in 
Fluid Particle Dynamics M and in the Smoothed Profile 
method [32|. 

Once one settles on solvers for the electrostatic and hy¬ 
drodynamic equations, the next step is to choose a sim¬ 
ulation model for the large particle. In the case of poly¬ 
electrolytes, a common choice is the bead-spring model, 
which represents a polyelectrolyte as a chain of consec¬ 
utive charged beads connected by springs. Large rigid 
particles, such as micrometer-sized colloids, can be in¬ 
troduced through appropriate boundary conditions, e.g., 
no-slip boundaries for the hydrodynamics and hard im¬ 
penetrable boundaries for solvent and ion particles. An 
alternative approach suitable for smaller colloids is the 
so-called “raspberry model”, which represents the col¬ 
loid by a shell made of of surface beads. The shape of 
the shell is maintained either by springs [~33| or by fixing 
the bead position with respect to the colloid center [34j . 

The computational expenses for the electrostatic and 
hydrodynamic computations are not equivalent. At high 
salt concentrations, a comparison fdhl has shown that the 
computational time is dominated by the costs of treat¬ 
ing the charges. Unfortunately, as mentioned above, im¬ 
plicit models for small ions are not suitable for dynamic 
studies. We recently proposed an efficient algorithm that 
overcomes the bottleneck caused by the explicit charges 
[36j |. It makes a compromise between computing effi¬ 
ciency and taking full consideration of correlations at all 
scales. The evolution of the ionic concentration is com¬ 
puted using Brownian pseudo particles [3^ , which is very 
fast because pseudo particles have no direct pair inter- 
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actions. In this approach, one chooses a coarse-grained 
length scale. Short-range correlations between small ions 
on length scales shorter than the coarse-grained length 
are neglected, but long-range charge-charge correlations 
can be retained. The full method consists of the pseudo¬ 
ion solver for electrolytes, DPD for the fluid and a P3M 
Coulomb solver. Simulations of electro-osmosis showed 
that the computer time required for the electrostatic cal¬ 
culations can be reduced to about half the time required 
for treating the fluid. Moreover, since the number of 
pseudo particles can be chosen independently of the num¬ 
ber of ions, it is independent of the ion concentration. 
Therefore, the proposed method is particularly suited for 
electrolyte solutions at high salt concentrations. 

Another issue in electrostatic simulations is the dielec¬ 
tric constant, which is often assumed to be homogeneous 
throughout the simulation box. In reality, the dielec¬ 
tric permittivity in the water and in the nanoparticles 
differ significantly. The difficulty to include dielectric 
contrast in molecular dynamics lies in the calculation of 
induced polarization charge at the interfaces, which has 
to be done self-consistently in every time-step. Recently, 
a number of authors have addressed this problem i38H40l ] . 
For example, Barros et al. [Ic} proposed to use a General¬ 
ized Minimum Residual method (4j| to calculate the dis¬ 
tribution of surface bound charges on nanoparticles and 
applied this approach to study the self-assembly of bi¬ 
nary colloids. They observed an unexpected string struc¬ 
ture formation due to the dielectric effect [42j . Boundary 
methods such as that sketched above also open up the 
possibility to study phenomena associated with induced- 
charge electrokinetics [43J. However, they cannot be used 
to simulate media with a smoothly varying dielectric con¬ 
stant e(r). Such systems can be treated using the re¬ 
cently proposed Maxwell Equations Molecular Dynamics 
(MEMD) method [2C],HI}, which is based on the Maggs 
approach [l9[ and solves the Maxwell equations by prop¬ 
agating a set of virtual auxiliary fields on a local scale. A 
recent review on methods to deal with dielectric contrasts 
can be found in Ref. [22| . 


III. PARTICLES IN CONSTANT ELECTRIC 
FIELDS 

Let us now consider an object (molecule or colloid) 
with a characteristic size of a suspended in an electrolyte 
solution. The object can acquire surface charges by sev¬ 
eral means, either by the dissociation of protons or by the 
selective adsorption of ions from the aqueous solution. 
Due to the electrostatic interaction, it is surrounded by 
a so-called “electric double layer” of oppositely charged 
ions. Some of these may bind strongly to the surface 
(Stern layer), and their main effect is to reduce the effec¬ 
tive surface charge. The others form the “diffuse layer”, 
which is characterized by a constant turnover of weakly 
bound, mobile ions. The thickness of the diffuse layer re¬ 
sults from the competition of osmotic pressure and elec¬ 


trostatic interactions and is typically of the order of the 
Debye screening length k . Both k - 1 and the particle 
extension, a, represent important characteristic lengths 
of the system. When an external electric field E is ap¬ 
plied, a positively charged particle starts to move in the 
direction of the electric field, and the ions in the diffuse 
layer experience a force in the opposite direction. Thus 
the particle experiences a friction force exerted by the 
electrolyte fluid which prevents its movement. In the 
stationary case, the electric driving force and the friction 
force balance each other and the object moves with a con¬ 
stant velocity v. The ratio between the terminal velocity 
v and the applied electric field E is called electrophoretic 
mobility /./, 

v = /rE. (1) 

In general, the electrophoretic mobility is a tensor, but 
we shall here restrict the discussion to isotropic objects, 
in which case the mobility is a scalar quantity. 

Before we proceed to presenting simulations, we briefly 
recapitulate the main contributions to the friction force 
involved in electrophoresis [44] 

1. The viscous force F v j s exerted by the fluid. The 
magnitude of this force is given by the Stokes fric¬ 
tion F v - ls = —Girr/av, where r] is the shear viscosity 
of the fluid. 

2. The electrophoretic retardation force F re t due to 
the movement of the counterion cloud. Since the 
counterions have the opposite charge, they move in 
the opposite direction of the central object. Ideally, 
the counterions compensate the charged object and 
the whole system is charge neutral. Therefore, the 
electrostatic force on the counterions exactly can¬ 
cels the driving force on the object, but to which 
extent this force is transferred to the central object 
is complicated. 

3. The polarization force F po i due to the distortion of 
the counterion cloud. As an example, let us con¬ 
sider a spherical colloid. The charge centers of the 
colloid and its surrounding counterions coincide at 
zero external field. If an external field is applied, 
the two centers are displaced slightly as the coun¬ 
terion cloud is distorted, which induces a dipole 
moment. The central object experiences an extra 
electric force due to this distortion. 


A. Electrophoresis of Colloidal Particles 

In general, the electrophoretic mobility of colloids de¬ 
pends on several parameters, e.g., the colloidal size a, the 
Debye screening length k , and the surface potential at 
the plane of shear, termed the ((-potential. In simulation 
studies, it is sometimes easier to prescribe the surface 
charge density er. Unfortunately, the relation between 
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a and £ is not simple j45| ; analytic formulas only exist 
for simple geometries (planes and cylinders). The elec¬ 
trophoretic mobility has a simple form only for special 
cases such as weakly charged colloids ((/fcsT C 1). In 
this limit, the counterion distribution can be treated us¬ 
ing the Debye-Hiickel approximation. Two well-known 
results are the Hiickel and Smoluchowski formulae: If 
the electric double layer is thick [kg <C 1), one can ne¬ 
glect the retardation and polarization forces, and write 
the electrophoretic mobility as J4(| 

M=|-C («a<l), (2) 

3 ry 

where e is the permittivity of the solution. This formula 
is derived from the balance between the electric driving 
force and Stokes friction. 

In the opposite limit of thin electric double layers, one 
obtains the Smoluchowski formula [13] 

H = -C («a > 1). (3) 

V 

In this case, the retardation force is accounted for but the 
polarization force is neglected. The situation is different 
if the surface charge is kept constant. In this case, the 
surface potential vanishes for large Ka, resulting a zero 
mobility. 

For intermediate values of Ka, Henry derived an ana¬ 
lytic formula of the form [48] 

H = f{na)~ C (4) 

with a scaling function /(kcs) that interpolates between 
the Hiickel and Smoluchowski result in their correspond¬ 
ing limits: f(na) —> 1 for Ka <C 1 and /(/ta) —>• 2/3 for 
kq^> 1. 

The assumptions and approximations made in deriving 
these simple formulae are: 

• The distortion of the counterion cloud is assumed 
to be negligible, and the contribution from the po¬ 
larization force Fp 0 i is omitted. 

• The distribution of salt ions and counterions is 
treated at the mean-field level using a continuum 
approximation. In equilibrium, this corresponds 
to the Poisson-Boltzmann approximation. Short- 
range correlations between charged species are ne¬ 
glected. 

• The colloidal particle is assumed to be weakly 
charged (e^/ksT < 1). This assumption justifies 
the use of the Debye-Hiickel approximation, and 
the simplification of the Poisson-Boltzmann equa¬ 
tion to a linear equation facilitates the derivation 
of analytic expressions. In this limit, the counte¬ 
rion distribution near the charged surface decays 
exponentially with the characteristic length k -1 . 


We shall now discuss how theoretical considerations 
and numerical calculations in the last few decades have 
helped to probe phenomena which are beyond those ap¬ 
proximations. 


1. Methods Based on the Electrokinetic Equations 

A general framework to study the dynamic phe¬ 
nomenon in electrolyte solutions is provided by the elec¬ 
trokinetic equations [lj. This involves a change in the 
point of view: Instead of calculating the friction force on 
the charged colloids, one writes down explicitly a set of 
partial differential equations which governs the dynam¬ 
ics of electrolyte solutions at the level of a continuum 
theory. The electrokinetic equations have three basic 
constituents: The Poisson equation for the electrostatic 
potential, the Navier-Stokes equations characterizing the 
fluid flow, and the Nernst-Planck equation, which is basi¬ 
cally a convection-diffusion equation describing the time 
evolution of ionic concentration profiles. Since thermal 
fluctuations and the discrete character of charges are ne¬ 
glected in this continuum approach, it corresponds to a 
mean-field approximation and reproduces the Poisson- 
Boltzmann results at equilibrium. All components of the 
friction force are implicitly included in this approach. 

In an early seminal paper [4f| , O’Brien and White com¬ 
puted the electrophoretic mobility of a single sphere in 
an infinite domain by numerically solving the electroki¬ 
netic equations for weak external fields (linear regime). 
Their results differed qualitatively from the previous pre¬ 
dictions in the case of the thin electric double layer: In¬ 
stead of being a monotonically increasing function of the 
zeta-potential £, as predicted by Eq. m , the mobility 
has a maximum for values of Ka > 3. This is shown in 
Fig-11! right. The non-monotonic behaviour results from 
the fact that different competing factors contributing to 
electrophoresis scale differently with £: The electrostatic 
driving force is proportional to £, while the friction force 
associated with the distorted counterion cloud is propor¬ 
tional to £ 2 . 

Based on the same electrokinetic model, Schmitz and 
Dunweg f50| recently developed a lattice-based approach 
to solve the linearized electrokinetic equations. Their 
strategy is divide-and-conquer: They divided the original 
equations into parts, and developed different numerical 
solvers for each equation. They then combined different 
solvers to compute the original equations using an iter¬ 
ation procedure. Their method is different to O’Brien 
and White in two aspects: One is the usage of periodic 
boundary condition, which allows to study the effect of 
colloid concentration. The second difference is that no 
assumption is made about the colloid shape, and one can 
in principle study colloidal particle with irregular shape. 

It is also possible to start with the full electrokinetic 
equations instead of the linearized version. The three 
equations belong to different categories of partial differ¬ 
ential equations, and there exist efficient solvers for each 
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FIG. 1: Reduced mobility E = ( 6 ^ 77 e/ek B T)^i as a 
function of the reduced zeta potential y = (e/k B T)£. 
Left: The mobility is a monotonic function for thick 
electric double layers na < 3. Right: A mobility 
maximum appears for thin electric double layers, 
na > 3. Reprinted from Ref. [49| with permission from 
The Royal Society of Chemistry. 


equation separately. Based on the Lattice Boltzmann 
method, Giupponi and Pagonabarraga studied the elec¬ 
trophoretic mobility of charged colloids 5l|. The Nernst- 
Planck equation was solved using a discrete method [52}. 
The results match those of O’Brien and White well in 
the case of small zeta-potentials, but deviate when the 
zeta-potential is large. In this regime, the method used 
by O’Brien and White suffers from numerical problems 
[49(, which may explain the discrepancy. In particular, 
Giupponi and Pagonabarraga found that a mobility max¬ 
imum exists for all salt concentrations. They further in¬ 
vestigated the effect of diffusivity of the small ions. 


Kim et al. took a different approach to the fluid by 
solving the Navier-Stokes equations directly 53]. The 
difficulty lies in the proper and efficient treatment of 
the moving boundary of the colloids. They devised 
the Smoothed Profile method, which replaces the shape 
boundary by a smooth interface with finite thickness 
f32j . Their results for the electrophoretic mobility showed 
good agreement with O’Brien and White for thick electric 
double layer (na = 0.5). The efficiency of their numerical 
method permits the simulation of many colloids in one 
simulation box. They further investigated the mobility 
dependence on the colloid concentration and compared 
their results with the theoretical prediction of Ohshima 
54]. At low colloid concentration, the agreement be¬ 
tween the simulation and theory is good. Deviations 
become noticeable when the electric double layers from 
different colloids overlap. 


2. Particle-based, Simulations 

All studies based on the electrokinetic equations ne¬ 
glect short-range correlations, as the small ions are 
treated in terms of ionic concentration fields. To remove 
this approximation, one must simulate the small ions as 
particles with excluded volume interactions that carry 
discrete charges. Combined with mesoscopic methods for 
the fluid, one obtains a simulation scheme which accounts 
for correlations between charged species. One drawback 
is an increase in computational time which makes the 
simulation of large colloids very difficult. 

Lobaskin et al. studied the electrophoresis of charged 
colloids in electrolytes containing only counterions or 
with very low salt content [55l . [56}. In their approach, 
the fluid is simulated using the LB method, and the col¬ 
loidal particle is modeled using the raspberry model [13 ■ 
Using a dimensional analysis, they demonstrated that the 
reduced electrophoretic mobility, defined by fj , le d = /V/^h 
(^h is the Hiickel result, Eq. [2j, depends only on two di¬ 
mensionless parameters Z e sl B /a and na. Here Z e s is 
the effective charge of the colloid, which differs from the 
bare colloid charge Z by the amount of strongly adsorbed 
counterions, Ib = e 2 / 4ir£k B T is the Bjerrum length, and 
the inverse screening length n is defined as 

K 2 = A-Kl B {ci + c s ), (5) 

where Cj = Z/U is the counterion concentration and c s 
is the ionic concentration due to the salt ions. Based 
on the dimensional analysis and Eq. © , one can map 
the mobility of colloids in salt-free fluids (the number of 
counterions is set by the colloid concentration) to that 
in electrolytes containing a small amount of salt at the 
same value of na. This correspondence is confirmed by 
both the simulations and experiments [55} . 

A similar model was used by Chatterji and Horbach in 
a series of studies of highly charged colloids [H, H2, [58} . 
They varied the colloid charge density instead of the 
zeta potential because the former can be controlled more 
easily in simulations. Upon increasing the charge den¬ 
sity, the electrophoretic mobility was found to initially 
increase, then reach a maximum, and decrease again. 
These simulation results are in accordance with the pre¬ 
dictions for the thin electric double layer case (Fig. [I] 
right). 

Chatterji and Horbach also considered systems con¬ 
taining divalent counterions, in which they found the 
electrolyte mobility to be reversed at high charge den¬ 
sity, indicating an overcompensation of the surface charge 
by the multivalent counterions. This overcharging phe¬ 
nomenon can be explained by ion correlations [591 ] and 
has also been observed experimentally. In experiments, 
however, the effect is sometimes larger than predicted by 
numerical simulations, suggesting that it might be en¬ 
forced by ion specific attractive forces [gQ]. Recent simu¬ 
lations by Raafatnia et al. have shown that overcharging 
may even occur in electrolytes containing only monova¬ 
lent ions if the colloids are coated by a suitable organic 
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layer j6li, [62]. 

As mentioned above, treating colloids with sizes much 
larger than the ion size is difficult with explicit ion mod¬ 
els. Raafatnia et al. have developed a hybrid approach to 
studying the electrophoresis of large colloids [63|. They 
first measured the zeta-potential of a flat surface us¬ 
ing simulations with explicit microions. The measured 
zeta-potential was then used as an input in the electroki- 
netic equations to calculate the electrophoretic mobility. 
When compared with experiments, the scheme worked 
well for monovalent and divalent salt solutions. In the 
case of trivalent salt, they needed to introduce an at¬ 
tractive interaction in order to reproduce the mobility 
reversal which is observed experimentally. 


B. Polyelectrolytes 

We turn to discussing polyelectrolyte electrophoresis. 
Compared to colloid electrophoresis, there are several dif¬ 
ferences: 

• Colloidal particles are usually at least one order 
of magnitude larger than the small molecules (sol¬ 
vents and salt ions). The separation of length scales 
makes it possible to describe large colloids in terms 
of boundary conditions and motivates the use of 
continuum approaches. In the case of polyelec¬ 
trolytes, the size of the monomer units is compa¬ 
rable to that of small molecules and the applica¬ 
tion of boundary methods becomes questionable. 
Even though the zeta potential is sometimes used 
to parameterize experimental results, its physical 
meaning is not always obvious. In simulations, an 
explicit treatment of small ions is often more ap¬ 
propriate than the use of continuum theories. 

• Whereas colloidal particles are rigid objects, poly¬ 
electrolyte chains can assume many conformations. 
The conformational dynamics is strongly influenced 
by the electrostatic and hydrodynamic interactions 
between monomers. Therefore, a full treatment of 
both interactions is necessary. Implicit methods 
have been designed to circumvent this requirement 
for stationary situations [64ll65j| , but in general, one 
has to be careful when dealing with non-stationary 
states. The conformational flexibility complicates 
studies of electrophoretic mobility. 

• If polyelectrolytes are highly charged, the strong 
electrostatic interaction attracts counterions to 
the proximity of the chain backbone. This phe¬ 
nomenon, called Manning condensation [6f|, also 
influences the electrophoretic mobility. 

• In the case of colloidal particles, the size and the 
surface charge density can be adjusted indepen¬ 
dently. For polyelectrolytes, this is usually not 
the case. The control parameters are the chain 


length N and the charge fraction. We shall re¬ 
strict ourselves to strongly charged chain, where 
each monomer carries a unit charge. The size of a 
polyelectrolyte chain is characterized by the radius 
of gyration R gi which depends on the chain length. 
Therefore, the electrophoretic mobility is usually 
represented as a function of the chain length. 


1. Free Solution Electrophoresis 

When the size of polyelectrolyte chain is much smaller 
than the Debye screening length, nR g -C 1, one may 
invoke the Hiickel picture for colloidal particles and con¬ 
sider only the electric driving force and the viscous fric¬ 
tion. The driving force of the external field is propor¬ 
tional to the chain length, while the viscous force from the 
fluid also increases with the chain length, but to a smaller 
extent due to the hydrodynamic interactions. Monomers 
experience hydrodynamic drag and shield each other 
from the exposure to the outer fluid, which effectively 
reduces the total Stokes friction. Therefore, for short 
chains, the mobility increases with the chain length. 

This behaviour was indeed observed in simulations and 
experiments. Grass and Holm modeled the polyelec¬ 
trolyte as a chain of charged beads connected by finitely 
extensible nonlinear elastic bonds [67]. The fluid was 
simulated using the LB method, and the Coulomb inter¬ 
action was calculated using the P3M method. They com¬ 
puted the electrophoretic mobility from equilibrium sim¬ 
ulations using the corresponding Green-Kubo relation. 
Later they also directly measured the terminal velocity 
in an external electric field [H, [6^], and found that the 
results agree in the weak-field regime. Frank and Win¬ 
kler used a similar model for the polyelectrolyte chain, 
but simulated the fluid using MPCD (70j |. Both studies 
emphasize the importance of hydrodynamic interaction, 
as complementary Langevin simulations show a decrease 
in the mobility with increasing chain length. 

One peculiar observation in these simulations is the 
presence of a mobility maximum at chain lengths N ~ 10 
(see Fig. [2]). This was attributed to a reduction of the ef¬ 
fective charge of the polyelectrolyte chain, which reduces 
the electric driving force. As the chain becomes longer, 
the polyelectrolyte assumes a rod-like conformation due 
to the mutual repulsion of monomers, and Manning con¬ 
densation sets in. The accumulation of counterions on 
the backbone reduces the effective charge of the chain, 
thereby reducing the electric driving force. Frank and 
Winkler measured the number of condensed counterions 
in their simulation using a distance criterion [7Cj. They 
found that as the chain becomes longer, the ratio between 
the condensed counterion to the total charge of polyelec¬ 
trolyte chain approaches 1 — 1/C, where £ is the Manning 
parameter. Grass and Holm applied different estimators 
for the effective charge, with similar results [§8]. How¬ 
ever, the exact point at which the saturation is reached is 
still unclear, because the polyelectrolyte chain is flexible, 
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N 


FIG. 2: Reduced electrophoretic mobility fj. of a 
polyelectrolyte chain as a function of the chain length 
N. The mobility is scaled with the free-draining limit 
^fd- Reprinted from Ref. [68j] with permission from 
The Royal Society of Chemistry. 


whereas the theory of Manning condensation applies for 
rigid rods. 

When the chain length increases further, the elec¬ 
trophoretic mobility reaches a constant value. This so- 
called “free-draining” limit has been well studied in the 
literature 0, Q- The fact that the mobility does not 
depend on the length of the polymers prevents the sepa¬ 
ration of long polyelectrolyte chains by free solution elec¬ 
trophoresis. The physics behind the plateau in mobility 
is two-fold: On the one hand, the effective charge per 
monomer becomes a constant, once the Manning con¬ 
densation sets in. Therefore, the electric driving force 
is proportional to the chain length. On the other hand, 
the effective friction force opposing the electric force also 
scales linearly with the chain length once the polyelec¬ 
trolyte size becomes larger than the Debye screening 
length. This is because a long chain can be viewed as 
a chain of charge neutral blobs with the blob size being 
the Debye length. When applying an external electric 
held, the chain experiences a force, but the surrounding 
counterions experience a force in the opposite direction, 
such that the sum of all forces on charged particles in a 
blob is zero. Therefore, no how is induced, the hydrody¬ 
namic interactions associated with the electric force are 
screened [72], and the resulting electrophoretic mobility 
does not depend on the chain length. 


2. Effect of Adding Salt 

The hydrodynamic screening effect discussed above be¬ 
comes even more pronounced in the presence of salt. 
Salt plays a dual role in polyelectrolyte electrophore¬ 
sis, because it screens both the electrostatic and the hy¬ 
drodynamic interactions. The screened electrostatic po¬ 


tential decays exponentially with a characteristic length 
/■v 1 . The screened hydrodynamic how profile still fea¬ 
tures a long-range 1 /r 3 decay, which is due to the vecto¬ 
rial nature of the velocity held [73[. In free solution, the 
screening length for hydrodynamics At^ 1 coincides with 
the Debye screening length At -1 . Therefore, varying the 
salt concentration induces a complicated interplay be¬ 
tween the electrostatic and hydrodynamic interactions. 
In Ref. [69], Grass and Holm examined the effect of salt 
on the electrophoretic mobility of short polyelectrolyte 
chains. The increase in the salt concentration reduces 
the screening length /tj^ 1 and therefore reduces the ini¬ 
tial shielding of hydrodynamics. As a result, the elec¬ 
trophoretic mobility becomes almost length-independent 
at high salt concentration. Fischer et al. [74] found that 
the counterion mobility changes sign as a function of salt 
concentration. At low salt concentration, the counterions 
condense to the polyelectrolyte backbone and the hydro- 
dynamic drag forces the counterions to move together 
with the polyelectrolyte chain. At high salt concentra¬ 
tions, however, the screening of hydrodynamics decou¬ 
ples the motion of counterions and polyelectrolyte, and 
counterions move in the opposite direction of the poly¬ 
electrolyte. 


3. Effect of Confinement 

Another way of decoupling the screening lengths for 
electrostatics and hydrodynamics is to use confinement. 
The screening length for hydrodynamics is given by ei¬ 
ther the characteristic size of the confinement or the De¬ 
bye length, depending on which is smaller. Hickey et al. 
investigated the electrophoretic stretching of a polyelec- 
trolyte chain confined between parallel plates [751] . They 
found that the hydrodynamic interaction is screened by 
the confining walls for strongly confined chains. Their 
study focused on uncharged walls. In the case of charged 
wall, the electroosmotic flow induced by the counteri¬ 
ons from the charged surface can reverse the movement 
of the polyelectrolyte chain: a positively charged poly¬ 
electrolyte chain can move in the opposite direction of 
the external held if confined by negatively charged walls 

mm. 


IV. PARTICLES IN ALTERNATING ELECTRIC 
FIELDS 

Alternating electric fields (AC Helds) provide another 
attractive tool for manipulating particles. Compared to 
constant electric Helds (DC fields), one has several ad¬ 
vantages: 

• In homogeneous AC fields, the displacement of the 
charged particle averages to zero. This is in con¬ 
trast to the DC case, where particles may travel 
over long distances during the time of an experi- 








ment. Since charged particles are mostly station¬ 
ary, electrodes can be placed close to each other. 
Therefore, one can produce large value of the elec¬ 
tric field, which is in general difficult in the DC 
case. 

• Using AC fields, one can avoid the accumulation of 
charged species on electrodes. 

• Unlike DC fields, AC fields do not generate con¬ 
stant electro-osmotic flows. 

• Apart from the amplitude, one can also tune the 
frequency and the phase of AC fields. Since dynam¬ 
ical processes in the system can take place on dif¬ 
ferent time scales, a time-dependent perturbation 
can probe the dynamics on selective time scales. 

• Finally, AC fields can also be used to manipulate 
uncharged particles. This is because they can in¬ 
duce polarization charges both in charged and un¬ 
charged particles, and the polarized particles can 
then be further manipulated by applying electric 
field gradients. 


A. Colloidal Particles 


1. Maxwell-Wagner Theory and Electrokinetic Theory 

The calculation of the induced dipole moment of a sin¬ 
gle spherical particle immersed in the medium provided 
by an electrolyte solution is a classical problem in elec¬ 
trodynamics, known as Maxwell-Wagner mechanism of 
dielectric dispersion jH, [§6j. The complex polarizability 
a(ui) has the form 


a(uj) = 47T£ m a 3 /C(w) = 47re n 


^ + 2e* 


(V 


where fC(u>) is the Clausius-Mossotti factor. The sys¬ 
tem is characterized by the complex dielectric functions 
6p m = £p,m + Kp,m/iw, where e is the permittivity and 
K the conductivity, and the subscript (p, m) labels the 
particle or the medium. From the form of the complex 
dielectric constant, one notices that the permittivity con¬ 
tribution dominates at high frequency (w —>• oo), while 
the conducting properties are more important at low fre¬ 
quency (lo —> 0). The frequency separating these two 
limits is the inverse of the Maxwell-Wagner relaxation 
time 


7m w 


“h 2<S m 
K p + 2K m 


( 8 ) 


The dielectric response to AC fields can be character¬ 
ized by the polarizability a.{w), where w is the field fre¬ 
quency. In analogy to the electrophoretic mobility, the 
polarizability is defined as the ratio between the induced 
dipole moment and the external field, 

p = cn(w)E. (6) 

The dipole moment has contributions from the colloidal 
particle and its surrounding electric double layer. In gen¬ 
eral, the polarizability depends on both the frequency 
and the amplitude of the external field. In weak fields, 
however, linear response applies and the polarizability 
does not depend on the held strength. In strong electric 
Helds, nonlinear effects may become important. 

One important application of AC Helds is dielec¬ 
trophoresis [zl, [79j. In dielectrophoresis, the time- 
averaged force on a particle has the form 5ft{a(w)}V|E| 2 , 
where 3?{a(aj)} is the real part of the complex polariz¬ 
ability. Under the influence of a position-dependent AC 
field, the particle is driven along the direction of the field 
gradient [80 h 82||. The magnitude of the force depends 
on the polarizability of the particle, which permits sepa¬ 
ration of colloidal particles and biomacromolecules [83j |. 
In a spatially homogeneous setup, AC fields can be used 
to control the self-assembly of many particles [s3 |. The 
control is realized by tuning the induced dipole-dipole in¬ 
teraction between particles. The most important quan¬ 
tity in this context is the complex polarizability of the 
charged particle. Here we shall focus on studies of single 
colloidal particle and its polarizability under a homoge¬ 
neous field. 


For charged colloids, the electric double layer in¬ 
troduces another contribution to the dipole moment. 
O’Konski proposed a surface conductance term to ac¬ 
count for the effect of electric double layer {83], and the 
same mechanism was also extended to the study of el¬ 
lipsoidal colloids J88|. The Maxwell-Wagner-O’Konski 
theory has a simple analytic formulation which gives 
qualitatively correct predictions, but it suffers from one 
main drawback: Since the theory is entirely formulated 
in terms of macroscopic properties, such as the conduc¬ 
tivities and permittivities, the polarization charges are 
taken to be localized at the particle/nredium interface, 
and their spatial distribution is entirely neglected. This 
simplification is only justified for thin electric double lay¬ 
ers and in the high-frequency regime. 

Experiments revealed another dispersion in the low- 
frequency regime, which cannot be explained by the 
Maxwell-Wagner-O’Konski theory. This low-frequency 
dispersion, often called alpha-relaxation, results from the 
distortion of the electric double layer. To properly treat 
the electric double layer, one can resort to the electroki¬ 
netic equations. Dukhin and Shilov made the further 
assumption that the electric double layer is at local equi¬ 
librium with the surrounding bulk solution [89j] , which is 
only valid in the low-frequency regime. They derived an 
analytic theory for the thin electric double layer case, 
which correctly predicts the low-frequency dispersion. 
The theory has been extended to asymmetric electrolytes 
}9Cl - !93[ and to aspherical colloids 0. 

In situations that involve thick electric double layer 
and the whole frequency spectrum, one can solve the 
full electrokinetic equations using a variety of numerical 




9 


methods. DeLacey and White [95[ extended the method 
of O’Brien and White (HJ to AC fields and presented nu¬ 
merical solutions for the polarizability of a single spher¬ 
ical colloid [95]. The electrokinetic equations are lin¬ 
earized in term of the external field, thus the calculation 
assumes weak fields. The method was further refined by 
Mangelsdorf et al. [H, [97| and Zhou et al. Hill 

et al. developed an alternative numerical scheme which 
overcomes a numerical instability in the high-frequency 
regime .‘MIC an d ex tended the method to polymer-grafted 
particles [lOdl llOlj . The effect of the colloid concentra¬ 
tion can also be incorporated by performing the calcula¬ 
tions in a sphere, whose volume is eq ual to the inverse of 
the colloid number density 102-104j. Shih et al. applied 
the Smoothed Profile method [53 to study the response 
of charged colloids to AC fields and reproduced both 
the Maxwell-Wagner relaxation and the alpha-relaxation 

fioi noa. 

The numerical methods were also extended to rodlike 
particles. Zhao investigated the case of long parallel rod 
with two spherical caps on both ends [1071 ]. The calcu¬ 
lation showed that the transition frequency in the low- 
frequency dispersion is reduced when the rod becomes 
longer, and do not change once the rod length reaches a 
certain limit. Dhont and Kang developed theories based 
on electrokinetic equations to study two special cases: 
Very weakly charged rods where the dipole mom ent i s 
mainly induced by the obstacle of the solid colloids |108[ , 
and strongly charged rods where con dense d counterions 
fully compensate the bare rod charge jl09f . 


2. Simulations with Explicit Ions 


The present authors have used particle-based method 
to investigate the diel ectri c response of a charged 
nanometer-sized colloid jlldl4ll3| . The fluid dynamics 
is simulated using DPD, where small ions are included 
explicitly as charged beads with excluded-volume inter¬ 
actions. The large colloid is modeled using the raspberry 
model. We systematically investigated the complex po¬ 
larizability a{w} as a function of the frequency, shown in 
Fig. [3] 

Let us first discuss the situation in the low-frequency 
regime. In the absence of an external field, the system 
including the charged colloid and its surrounding electric 
double layer has spherical symmetry. When an external 
field is applied, the colloid (positively charged) moves in 
the direction of the electric field, while the counterions 
move in the opposite direction. This creates a dipole 
moment, which points in the same direction than the 
external field. It also distorts the electric double layer, 
compressing it at the front end and expanding it at the 
back end, and thus creates a concentration gradient from 
back to front. In addition, the colloid acts as a mov¬ 
ing obstacle for the salt ions in solution which are also 
driven by the electric field. Negatively charged ions ac¬ 
cumulate at the front end of the colloid, which combined 



FIG. 3: Real and imaginary part of the complex 
polarizability a(w) of a charged particle as a function of 
the frequency. The field strength is set in the linear 
region. The points with err-bars are simulation results. 

The solid lines give the prediction from the 
Maxwell-Wagner-O’Konski theory. The dashed lin es ar e 
numerical solutions to the electrokinetic equations (1011. 

Reprinted from Ref. [l 13j |. 


with the extra counterions results in a further increase of 
the salt concentration there, whereas the concentration 
at the back is further reduced. The concentration gra¬ 
dient induces a diffusive migration of the salt molecules 
from the front to the back. This concentration-induced 
effect reduces the dipole moment created by the external 
field. However, since it is a second-order effect caused 
by the field-induced dipole, the net dipole moment still 
points in the direction of the external field. This can be 
seen from Fig.[3l where the real part Re{a} has a positive 
value at low frequency. 

The diffusion of the salt over the distance of colloid 
diameter requires time, which can be estimated by 


t c ~ 


(2 a? 
A 


(9) 


where A is the diffusion constant of the small ions. If 
the field frequency increases beyond l/r c , the external 
field oscillates so fast that the diffusion process cannot 
follow. Therefore, the concentration-induced process is 
suppressed at frequencies / > l/r c , which effectively in¬ 
creases the dipole moment. This is the origin of the low- 
frequency dispersion and explains the slight increase of 
5i{a} in Fig. [3] 

When the frequency is further increased, the external 
field eventually oscillates so fast that the ion cloud cannot 
respond. Thus, both Re{a} and Im{a} drop to zero. The 
characteristic time corresponds to the time required for 
the small ions to diffuse over the distance of Debye length 


A 


K 


( 10 ) 
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The last equality applies for our simulation parameters 
(e m = £ p and K p = 0), which differs from Eq. © only 


by a prefactor. 

The simulation results can be compared with the theo¬ 
retical and numerical predictions. The Maxwell-Wagner- 
O’Konski theory gives the solid curves in Fig. [3] which 
are only in qualitative agreement with the simulation. 
The theory predicts roughly the correct transition fre¬ 
quency in the high-frequency regime, but misses the low- 
frequency dispersion. The numerical solution of the elec- 
trokinetic equations is shown with dashed lines. It is in 
quantitative agreement with the simulation results. 

We have also investigated the dielectric resp onse of 
uncharged colloids to external AC fields [ill] . The 
Maxwell-Wagner theory, Eq. 0, predicts a complex po¬ 
larizability, simply due to the fact that the medium is 
conducting and the particle is not. Physically, the par¬ 
ticle acts as an obstacle for the flow of charges induced 
by the applied field, such that charges of opposite sign 
accumulate at both ends. Since uncharged particles have 
no electric double layer, one does not expect additional 
contributions. Indeed, the simulation data were found to 
be in almost perfect agreement with the predictions of 
the Maxwell-Wagner theory. Fine details of the charge 
distributions can be understood quantitatively within the 
electrokinetic theory developed by Dhont and Kang [108 1. 


B. Polyelectrolytes 

Simulating the polyelectrolyte response to AC fields is 
a challenging task and only a few researchers have tack¬ 
led this problem. Most studies focused on the confor¬ 
mational change of a single polyelectrolyte chain when 
a strong field is applied. For polyelectrolyte chains that 
are sufficiently long, or in the case of high salt concen¬ 
tration, hydrodynamic interactions are strongly screened. 
Therefore, most simulation studies used a Langevin ther¬ 
mostat which neglects the hydrodynamics, but includes 
small ions explicitly in order to account for the Manning 
condensation. 

Liu et al. studied the unfolding and collapse of a flexi¬ 
ble polyelectrolyte under a sinusoidal electric field |ll4j . 
They first measured the critical field strength at which 
the chain undergoes a transition to the fully extended 
state. Their results confirmed a theoretical prediction by 

Netz urn nia : The critical field scales as E cr it ~ N x , 
where x = — 1/2 for collapsed chains and x = —3^/2 for 
non-collapsed chains with the Flory exponent v. They 
also estimated the relaxation time by measuring the cor¬ 
relation time of the end-to-end distance. For certain pa¬ 
rameters of the AC field, the polyelectrolyte chain exhibit 
a stretch-collapse cycle. The simulations indicated such 
a behaviour occurs under two conditions: (1) the field 
strength is larger than the critical strength and (2) the 
frequency is comparable or less than inverse of the chain 
relaxation time. Hsiao et al. examined a similar system, 
but in trivalent salt solutions under a square-wave elec¬ 


tric field [H3- They implemented a different estimator 
for the relaxation time. Since they use a square wave, 
the dipole moment exhibits an exponential decay after 
the electric field reverses its sign. They used the char¬ 
acteristic time for the exponential decay to estimate the 
relaxation time and reached the same conclusions than 
Liu et al. regarding the requirement for chain stretching. 


V. CONCLUSIONS 

Strategies to control particles with electric fields have 
attracted considerable attention in recent years. Since 
the particles respond to the external electric fields on rel¬ 
atively short time scales and in an often fully reversible 
way, using electric fields provides an attractive approach 
to control the position of individual particle or the struc¬ 
ture of particle assemblies. In this review, we have dis¬ 
cussed recent simulation studies of charged colloids and 
polyelectrolyte chains under external electric field. These 
studies have identified two important length scales: the 
Debye screening length ac _ 1 and the size of the particles 
(the radius a of colloidal particles or the radius of gyra¬ 
tion R g of polyelectrolytes). In the case of AC fields, the 
frequency-dependent response is determined by the diffu¬ 
sion times (r mw and r c ) associated with those two length 
scales. In most situations, the delicate interplay between 
the electrostatic and the hydrodynamic interactions plays 
an important role in determining the dynamics. 

Although many studies have addressed the behaviour 
of particles in external electric fields, many open ques¬ 
tions still remain. For example, we have only discussed 
the behaviour of spherical particles with homogeneous 
surface properties. Removing this constraint opens up 
the whole new territory of anisotropic particles Ill8[ . 
Colloidal particles can have heterogeneous surface charge 
distributi ons [m, mg, or varying hydrodynamic slip 
12lUl23l . These novel types of colloidal particles and 
their response to electric fields still remain to be explored. 

We have focused on flexible polyelectrolyte chains, 
while most interesting polymers are semiflexible, for ex¬ 
ample, double-stranded DNA. Understanding how semi¬ 
flexible polyelectrolytes move under external fields is not 
only important from a scientific point of view, but can 
also he lp ad vancing molecular biology and medicine re¬ 
search [124 |125| . Regarding the conformational changes 
of polyelectrolyte chains under external fields, we have 
only touched one aspect of the story: the chain stretch¬ 
ing due to the external fields. Experiments have shown 
that it is possible to col l apse polyelectrolyte chains un¬ 
der DC or AC fields }l26i . il27j . The mechanism for this 
surprising behaviour is still under debate, and input from 
theories and simulations would certainly help to under¬ 
stand such phenomenon. 

Furthermore, colloidal particles and polyelectrolytes 
chain can be seen as two extremes of particles with 
varying rigidity. There exists a zoo of so-called soft- 
particles that interpolate between the solid colloids and 
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flexible polyelectrolytes mm- Notable examples are 
polymer-grafted colloids, star polyelectrolytes, and mi¬ 
crogels. Simulations can help to understand how to con¬ 
trol and manipulate those particle using electric fields. 

In this review, we have focused on the case of dilute 
suspensions. In a dense suspension or at low salt con¬ 
centration, the electric double layers of different parti¬ 
cles may overlap. This introduces another length scale, 
the average distance between charged particles. External 
fields may modify the effective interaction among parti¬ 
cles, which in turn influences the self-assembly process 
in dense suspensions. This opens new possibilities for 
manipulating particles, which may be used for directed 
self-assembly. Simulations may be useful to guide the 


design of experimental protocols for making novel inter¬ 
esting nonequilibrium structures. 
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